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We briefly discuss the phenomenological theory of dissipative fluid. We also present some numer- 
ical results for hydrodynamic evolution of QGP fluid with dissipation due to shear viscosity only. 
Its effect on particle production is also studied. 
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I. INTRODUCTION 

A large volume of experimental data from Au+Au col- 
lisions at RHIC are successfully analysed in an ideal fluid 
dynamic model [lj. However, experimental data do show 
deviation from ideal behavior. The ideal fluid descrip- 
tion works well in almost central Au+Au collisions near 
mid-rapidity at top RHIC energy, but gradually breaks 
down in more peripheral collisions, at forward rapidity, 
or at lower collision energies @, indicating the onset of 
dissipative effects. To describe such deviations from ideal 
fluid dynamics, quantitatively, requires the numerical im- 
plementation of dissipative relativistic fluid dynamics. 

Though the theories of dissipative hydrodynamics 
@j S 111 has been known for more than 30 years , signif- 
icant progress towards its numerical implementation has 
only been made very recently @, 0, H, @, U3| • In the fol- 
lowing, I will briefly review the phenomenological theory 
of dissipative hydrodynamics. Some numerical results 
comparing the ideal and 1st order dissipative hydrody- 
namics, obtained using the computer code AZHYDRO- 
KOLKATA [3 will be shown. 



II. PHENOMENOLOGICAL THEORY OF 
DISSIPATIVE HYDRODYNAMICS 



In relativistic fluid dynamics, one defines a time-like 
hydrodynamic 4- velocity, (normalised as u 2 = 1). One 
also define a projector, A^ = g^ v — u^u v ', orthogonal to 
the 4- velocity {A^ v u v = 0). In equilibrium, an unique 
4- velocity exists such that the particle density (n), 
energy density (e) and the entropy density (s) can be 
obtained from, 
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An equilibrium state is assumed to be fully specified 
by 5-parameters, (n, e, u M ) or equivalently by the thermal 
potential, a = li/T (/j, being the chemical potential) and 
inverse 4-temperature, (3^ — u^/T. Given a equation of 
state, s = s(e,n), pressure p can be obtained from the 
generalised thermodynamic relation, 
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Using the Gibbs-Duhem relation, d(pfl^) — Njf q da — 
T^df3\, following relations can be established on the 
equilibrium hyper-surface T, eq (a, /3 M ), 
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In this section, I briefly discuss the phenomenological 
theory of dissipative hydrodynamics. More detailed ex- 
position can be found in [5J. 

A simple fluid, in an arbitrary state, is fully speci- 
fied by primary variables: particle current (A M ), energy- 
momentum tensor (T^) and entropy current (S r * t ) and 
a number of additional (unknown) variables. Primary 
variables satisfies the conservation laws; 



= 0, 
d^ v = 0, 



and the 2nd law of thermodynamics, 



d^S" > 0. 
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In a non-equilibrium system, no 4-velocity can be 
found such that Eqs I2.4I2.5I2.6I remain valid. Tensor de- 
composition leads to additional terms, 



s» = 
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The new terms describe a net flow of charge V* = 
A» V N V , heat flow, W = (e + p)/nV fl + g" (where 
is the heat flow vector), and entropy flow $ M . II = 
p is the bulk viscous pressure and ir^ = 
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tensor. Hydrodynamic 4-velocity can be chosen to elimi- 
nate either (the Eckart frame, is parallel to particle 



- A 1 " 7 A^ T - \A» V A° T ]T aT is the shear stress 
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flow) or the heat flow (the Landau frame, u M is parallel 
to energy flow). In relativistic heavy ion collisions, Lan- 
dau's frame is more appropriate than the Eckart's frame. 
Dissipative flows are transverse to u M and additionally, 
shear stress tensor is traceless. Thus a non-equilibrium 
state require 1+3+5=9 additional quantities, the dissi- 
pative flows II, (or V^) and tt^ . In kinetic theory, 
and are 1st and 2nd moment of the distribu- 
tion function. Unless the function is known a-priori, two 
moments do not furnish enough information to enumer- 
ate the microscopic states required to determine and 
in an arbitrary non-equilibrium state, no relation exists 
between, N v , T^ v and S^. Only in a state, close to a 
equilibrium one, such a relation can be established. As- 
suming that the equilibrium relation Eq |2.8l remains valid 
in a "near equilibrium state" also, the entropy current 
can be generalised as, 



S" = S£ + dS" = - aN» + /3 x T Xf * + Q» (2.13) 



where is an undetermined quantity in 2nd order in 
deviations, SN^ = A" 1 - Ng q and ST^ = T» v - Tj£. 
Detail form of Q M is constrained by the 2nd law d^S^ > 
0. With the help of conservation laws and Gibbs-Duhcm 
relation, entropy production rate can be written as, 
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As before, one can cast the entropy production rate 
(Td^S^) in the form of Eq |2.15l Neglecting the terms 
involving dissipative flows with gradients of equilib- 
rium thermodynamic quantities (both are assumed to be 
small) and demanding that a linear relation exists be- 
tween the dissipative flows and thermodynamic forces, 
following relaxation equations for the dissipative flows 
can be obtained, 
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where D — u^d^is the convective time derivative. Unlike 
in the 1st order theories, in 2nd order theories, dynami- 
cal equations controls the dissipative flows. Even if ther- 
modynamic forces vanish, dissipative flows donot vanish 
instantly. 



-SN^d^a + 5T^d^ u + d^Q* (2.14) 



Choice of leads to 1st order or 2nd order theories of 
dissipative hydrodynamics. In 1st order theories = 0, 
entropy current contains terms upto 1st order in devia- 
tions, 5N^ and 5T^ . Entropy production rate can be 
written as, 



TduS 1 * = HX - q^Xu + n^X 



(2.15) 



X^ = — u v d v u^ and X» v = 



where, X = — V.u 
S/<» u v> . 

The 2nd law, c^S 7 ' > can be satisfied by postulating 
a linear relation between the dissipative flows and and 
thermodynamic forces, 
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where £, A and r\ are the positive transport coefficients, 
bulk viscosity, heat conductivity and shear viscosity. 

In 1st order theories, causality is violated. If, in a given 
fluid cell, at a certain time, thermodynamic forces van- 
ish, corresponding dissipative fluxes also vanish instantly. 
This is corrected in 2nd order theories Q where entropy 
current contain terms upto 2nd order in the deviations, 
/ 0. The most general containing terms upto 2nd 
order in deviations can be written as, 
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FIG. 1: Contour plots of energy density at (proper) time 
r=5.2 fm. 

Before we proceed further, it may be mentioned that 
the parameters, a and j3\ are not connected to the actual 
state (N^,T^). The pressure p in Eq |2~T31 is also not 
the "actual" thermodynamics pressure, i.e. not the work 
done in an isentropic expansion. Chemical potential a 
and 4-inverse temperature (3\ has meaning only for the 
equilibrium state. Their meaning need not be extended 
to non-equilibrium states also. However, it is possible 
to define a fictitious "local equilibrium" state, point by 
point, such that pressure p in Eq |2.13l can be identified 
with the thermodynamic pressure, at least upto 1st order. 
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In 2nd order in deviations, such an identification is not 
possible. 
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FIG. 2: Contour plots of temperature at y=0 in x — r plane. 



III. VISCOUS HYDRODYNAMICS FOR QGP IN 
2+1 DIMENSIONS 

Numerically, causal dissipative hydrodynamics is a 
challenging problem. One needs to solve simultaneously 
14 partial differential equations(5 conservation equations 
and 9 additional equations for the dissipative flows:, II, 

and ■k^ v ). Recently, at the Cyclotron Centre, Kolkata, 
we have developed a code for solving causal dissipative 
hydrodynamics with dissipation due to shear viscosity 
only, in 1+1 dimension [8[ (assuming boost-invariance 
and cylindrical symmetry) and presently extending the 
code in 2+1 dimensions (with boost-invariance only). We 
have completed the coding for the 1st order theory. Here, 
we present some numerical results. 

Assuming boost-invariance, we have solved 1st order 
viscous hydrodynamics for initial state QGP in 2+1 di- 
mension, in (T,x,y,i] s ) coordinate. We restrict ourselves 
to central rapidity region, where the QGP fluid is essen- 
tially baryon free and to keep the calculations simple, 
we consider the most important dissipative term, the 
shear viscosity and neglect the other dissipative terms, 
e.g. heat conduction, bulk viscosity. For viscosity, we 
have used two values, the ADS/CFT motivated value 
?7/s=0.08, and the perturbative estimate ry/s=0.135 . De- 
tails of the equations solved can be found in [HI. We just 
mention that we have used the equation of state, EOS-Q, 
developed in ref.[l|, with 1st order phase transition with 
critical temperature 7^=164 MeV. As mentioned earlier, 
1st order dissipative hydrodynamics violate causality and 
can lead to unphysical effects like early reheating. But 
for QGP fluid, which is close to an ideal fluid, such effects 



can be minimised with appropriate initial conditions, and 
we did not find any evidence of early reheating. 

A. Evolution of the viscous fluid 

In the following we will show the results obtain in 
Au+Au collision at impact parameter b = 6.8 fm, which 
approximately corresponds to 16-24% centrality Au+Au 
collisions. With the same initial conditions, we have 
solved the energy-momentum conservation equations for 
ideal fluid and viscous fluid. In Fig[TJ we have shown the 
constant energy density contour plot in x-y plane, after 
an evolution of 5 fm. The black lines are for ideal fluid 
evolution. The red and blue lines are for viscous fluid 
with ?7/s=0.08 and 0.135 respectively. In Fig[21 contour 
plot of temperature in x-r plane is shown. It is evident 
that compared to ideal fluid, viscous fluid cools slowly. 
Transverse expansion is also enhanced in a viscous fluid. 
It is not shown, but fluid velocity grow faster in viscous 
flow. 
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FIG. 3: pr spectra of ir~ for ideal and viscous evolution. 

B. Particle spectra 

Viscosity influences the particle production by (i) ex- 
tending the freeze-out surface and (ii) by introducing a 
correction to the equilibrium distribution function (see 
[Irj||). In Fig'121 transverse momentum distribution of 
7T _ from ideal and viscous hydrodynamics are compared. 
Freeze-out temperature is Tp=0.158 GeV. Pion produc- 
tion is increased in viscous dynamics. We also note that 
effect of viscosity is more prominent at large px than at 
low px- Pt spectra of pions are flattened with viscosity. 

Effect of viscosity is also prominent on elliptic flow 
(Figf?]). In ideal dynamics, elliptic flow continues to in- 
crease with px- In viscous dynamics, on the otherhand, 
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FIG. 4: pr dependence of elliptic flow for n~ . 
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FIG. 5: pr spectra of n~ from ideal fluid with initial entropy 
density 110 fm~ 3 is compared with viscous dynamics with 
different initial entropy density. 

elliptic flow tends to saturate. The result is very en- 
couraging, as experimentally also elliptic flow tends to 
saturate at large pt- 



Lastly, in FigJSl we have shown a comparison of pt 
spectra obtained in ideal hydrodynamics with initial en- 
tropy density Si n i=110 fm~ 3 with px spectra obtained 
in viscous hydrodynamics obtained with initial entropy 
density, 60,80 and 110. pt spectra from viscous fluid ini- 
tialised with Si n i =60-80 compare well with the spectra 
from ideal fluid, initialised at higher entropy density. To 
produce the same spectra, viscous fluid require 20-30% 
less initial temperature. 

IV. SUMMARY AND CONCLUSIONS 

We have briefly reviewed the phenomenological theory 
of dissipative hydrodynamics and presented some numer- 
ical results from 1st order dissipative hydrodynamic of 
QGP fluid with only shear viscosity. 1st order theories 
suffer from the problem of causality, signal can travel 
faster than light. Unphysical effects like reheating of the 
fluid, early in the evolution, can occur. In this model 
study, we have considered two values of viscosity, the 
ADS/CFT motivated value, rj/s «0.08 and perturba- 
tively estimated viscosity, r]/s «0.135. We did not find 
any indication of unphysical reheating. Explicit simula- 
tion of ideal and viscous fluids confirms that energy den- 
sity of a viscous fluid, evolve slowly than its ideal counter- 
part. The fluid velocities on the other hand evolve faster 
in viscous dynamics than in ideal dynamics. Transverse 
expansion is also more in viscous dynamics. 

We have also studied the effect of viscosity on parti- 
cle production. Viscosity generates entropy leading to 
enhanced particle production. Particle production is in- 
creased due to (i) extended freeze-out surface and (ii) 
non-equilibrium correction to equilibrium distribution 
function. With ADS / CFT (perturbative) estimate of vis- 
cosity, at pr=3 GeV, pion production is increased by a 
factor 3 (5) . Increase is even more at large pr- While 
viscosity enhances particle production, it reduces the el- 
liptic flow. At p T =i GeV, for ADS/CFT(perturbative) 
estimate of viscosity, elliptic flow is reduced by a factor 
of 2(3). We also find that at large pr elliptic flow tends 
to saturate. 

To conclude, present study shows viscosity, even if 
small, can be very important in analysis of RHIC Au+Au 
collisions. Currently accepted initial temperature of hot 
dense matter produced in RHIC Au+Au collisions, ob- 
tained from ideal fluid analysis can be changed by 20% 
or more with dissipative dynamics. 
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